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Abstract. The special nature of gradient Yamabe soliton equation which was first observed by 
Cao-Sun-Zhang CSZ^ shows that a complete gradient Yamabe soliton with non-constant potential 
function is either defined on the Euclidean space with rotational symmetry, or on the warped 
product of the real line with a manifold of constant scalar curvature. In this paper we consider 
jy^ the classification in the latter case. We show that a complete gradient steady Yamabe soliton on 

warped product is necessarily isometric to the Riemannian product. In the shrinking case, we show 
that there is a continuous family of complete gradient Yamabe shrinkers on warped products which 
are not isometric to the Riemannian product in dimension three and higher. 
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^ 1. Introduction 

a 

Geometric flows are important tools to understand the topological and geometric structures 
in Riemannian geometry. A special class of solutions on which the metric evolves by dilations 
and diffeomorphisms plays an important role in the study of the singularities of the flows as they 
appear as possible singularity models. They are often called soliton solutions. In the case when the 
diffeomorphisms are generated by a gradient vector field, we call such soliton a gradient soliton. 
In this paper we are interested in the gradient soliton solutions to the Yamabe flow. This flow has 
been very well-understood in the compact case, see the very recent survey |Brej by S. Brendle and 
the references therein. For the non-compact case, in |DS1| P. Daskalopoulos and N. Sesum showed 
that the solutions to the Yamabe flow from some complete metric develop a finite time singularity 
and the metric converges to a soliton solution after re-scaling. 

A complete Riemannian manifold (M,g) is called a gradient Yamabe soliton if there exists a 
smooth function / : M — )■ R and a constant p£E such that the Hessian of / satisfies the equation 

(1.1) Ressf = (R-p)g, 

where R is the scalar curvature of M. If p > 0, p < or p = 0, then (M, g) is called a Yamabe 
shrinker, Yamabe expander or Yamabe steady soliton respectively. In [DS2J Daskalopoulos and 
Sesum initiated the investigation of gradient Yamabe solitons and showed that all complete locally 
conformally flat gradient Yamabe solitons with positive sectional curvature are rotationally sym- 
metric. This result is inspired by the work of the classification of locally conformally flat Ricci 
solitons, especially in |CaChj . 



The equation (1.1) where the right hand side is a smooth function, not necessarily the one given 



by the scalar curvature, was appeared in 1925 on the study of Einstein metrics that are conformal 
to each other by H. Brinkmann, see |Bri] . Such equation was also considered by J. Cheeger and 
T. Colding in their work [C hCoj . By exploring the special nature of the Yamabe soliton equation 



i 



2 



CHENXU HE 



( 1.1 ) H.-D. Cao, X. Sun and Y. Zhang showed that every complete gradient Yamabe soliton admits 
a global warped product structure in [CSZ]. 

Theorem (Cao-Sun-Zhang). Let (M n ,g,f) be a complete gradient Yamabe soliton satisfying equa- 
tion fPp with non constant function f . Then |V/| is constant on regular level hypersurfaces of f, 
and either 

(1) / has a unique critical point, and (M,g) is rotationally symmetric and equal to the warped 
product 

([0,Oo),dr 2 ) X| V /| {§ n -\gcan) , 

where g can is the round metric on S n_1 , or 

(2) / has no critical point and (M,g) is the warped product 

(M,dr 2 ) x, v/ | (N n -\g), 

where (iV n-1 ,g) is a Riemannian manifold of constant scalar curvature, say R. Moreover, 
if the scalar curvature R of (M,g) is non-negative, then either R > 0, or R = R = and 
(M,g) isometric to the Riemannian product (M, dr 2 ) x (N n ,g). 

By the global warped product structure from the theorem above, the Yamabe soliton equation 
( 1.1 ) reduces to the following ordinary differential equation in cp = | V/|, see [CSZ| Equation (2.17)]: 

/ / \ 2 // 



(1.2) <p' + p = <p - i R-(n-l)(n-2) — - 2(n - 1 

\<Pj ' ¥ 

Note that the scalar curvature R of M is given by R = <p' + p. In the first case when / has a unique 
critical point or (p(0) = 0, the above equation (1.2) is equivalently to equation (1.3) in [DS2] by 
changing variables. If one chooses the round metric g can on S" _1 with radius one, i.e., <p'(0) = 1, 
then R = (n — l)(n — 2) and p appears as a parameter in the differential equation of ip. The 
arguments for their equation (1.3) in |DS2j which can also be derived from equation (1.2) show 
that there is a unique complete Yamabe soliton metric for every p G M. The asymptotic behavior 
of (p which determines the asymptotic geometries in various cases can also be derived from their 
arguments. 

In this paper we consider the classification of Yamabe solitons in the second case, i.e., the manifold 
is topologically a product as M = K x iV™ -1 with (N,g) being of constant scalar curvature R, and 
g = dr 2 + tp 2 g defines a complete metric on M where <p E C°° (K) is a positive function. We call 
a Yamabe soliton metric is a product soliton if (p is a positive constant, i.e., M is isometric to the 
Riemannian product M x N. Note that in |CSZj a gradient Yamabe soliton is called trivial if the 
potential function / is constant. Our first result shows that a gradient steady Yamabe soliton is 
necessarily a product soliton. 

Theorem 1.1. Any complete gradient steady Yamabe soliton on M n = R x^ iV n_1 is isometric to 
the Riemannian product with constant cp and (N, g) being of zero scalar curvature. 

It has the following corollary by using the results above in the case when / has a critical point. 

Corollary 1.2. Up to scaling, a non-trivial complete gradient steady Yamabe soliton is either the 
product soliton or the unique one on W 1 with rotational symmetry. 
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On the other hand, by studying the solutions to the differential equation in tp we obtained a 
family of examples for complete gradient Yamabe shrinkers which are not product solitons. More 
precisely we have 

Theorem 1.3. Let p > be a constant and n > 3. Suppose R > p^x| (( n — l)( n + 2)) n + 2 is a 
positive constant and let (iV n_1 ,g) be a Riemannian manifold with constant scalar curvature R. 
Then there is a unique complete gradient Yamabe soliton metric on the warped product M n = 
RXy N' 11 ^ 1 with constant p that is not isometric to the product metric. 

We briefly discuss the proofs of our results. Since the variable r does not appear explicitly in the 



equation (1.2), by introducing new variables the equation can be reduced to a first order nonlinear 
differential equation. It turns out that the new equation is a classical one, the Abel differential 
equation of the second kind, see for example [In]. This equation does not have an explicit solution 
in general. So we consider the planar dynamical system defined by this equation. Every trajectory 
a with positive tp of the dynamical system defines a Riemannian metric on M. The completeness 
of the metric requires that certain integral divergences in both directions along a. In the steady 
case we show that all trajectories fail the divergence condition either in one or both directions. 
This gives us Theorem In the shrinking case, for each dimension n > 3, the integral test of the 
completeness singles out a unique trajectory 7 of this system. When R and p satisfy the equality 
in Theorem |1.3| instead of the strict inequality, the function tp is positive except for one point on 
7. On the other hand when R and p satisfy the strict inequality, by using some estimates on the 
trajectory 7 we can bound 7 by some simpler curves in the phase plane which allows us to show 
that tp is positive. 

Remark 1.4. The equality case of R and p suggests that ((« — l)(n + 2))™+ 2 is the optimal 
constant, i.e., we expect that any complete gradient Yamabe shrinker on M is isometric to the 



Riemannian product with tp being a constant function if the inequality of R and p in Theorem 1.3 
does not hold. 



Remark 1.5. It was pointed out in [DSS] that steady and shrinking Yamabe solitons have non- 
negative scalar curvature. From the second case in the theorem by Cao-Sun-Zhang, we can assume 
that R > in the proofs of Theorem 1 1 . 1 1 and Theorem 1.3 



Remark 1.6. In the expanding case, the Yamabe soliton equation defines dynamical systems that are 
similar to those in the shrinking and steady cases. The phase-plane analysis of such systems shows 
that there are complete Yamabe expanders on the warped products with non-constant function tp. 
For example when n = 6, for negative values of R there are such examples. The full classification 
in the expanding case will appear elsewhere. 



The paper is organized as follows. In Section 2 by using new variables the equation (1.2) in tp is 
reduced to a first order nonlinear differential equation and then we define the dynamical system for 
this first order equation. In Section 3 we analyze the system for steady Yamabe solitons and prove 



Theorem |1.1[ In Section 4 we study the system for Yamabe shrinkers and prove Theorem 1.3 In 



these two sections, one will see that the systems are considerably simpler when n = 6. So we prove 
the results for n = 6 first and then for other dimensions. There are several figures prepared by 
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the computer algebra system Mathematica which illustrate the ideas of the proofs in this paper. 
However our arguments do not rely on any of these figures. 

Acknowledgment. The author would like to thank Huai-Dong Cao for bringing the problems 
to his attention and many useful suggestions and discussions. He also wants to thank Xiaofeng Sun 
for helpful conversations. 



2. Preliminaries 



In this section we derive the first order differential equations from the equation (1.2) in ip and 
define a planar dynamical system for this equation. Then the problem of finding a complete gradient 
Yamabe soliton is equivalent to the one of finding a trajectory satisfies certain restrictions, see 



Proposition 2.1 



Since the variable r does not appear in equation (1.2) explicitly we let y = ip'(r) and then we 
have 



¥"{r) = °j^y = y'{p)y. 



Let 



x = ip(r) > 0, y = ip'(r), 

and then equation ( |1.2[ ) is given by 

2(n — l)xyy'(x) = R — (n — l)(n — 2)y 2 — x 2 (y + p). 

It can be rewritten as 



yy'{x) 



n-2 



-y 



R 



px 



2x a 2{n-iy 2{n-l)x 2(n-l) 
The above equation has the form of the Abel differential equation of the second kind 

yy'{x) = ai(x)y 2 + a 2 {x)y + a 3 (x) 

with 



Tl — 2 X 

ai{x) = — , a 2 (x) = -— — , a 3 (x) 



R 



px 



2x ' 2{n-l) w 2{n-l)x 2{n - 1) 

Using the standard process we write the equation in the canonical form. Let 



then we have 
with 

Fi(x) 



a 2 (x 



y = E(x)w, where E(x) = exp a\{x)dx 

ww'(x) = F\{x)w + Fq(x) 

a 3 (x) R 



_ n-2 
X 2 , 



1 



E{x) 2(n - 1 
By introducing the new variable 



-£2, F (x) 



Fi(x)dx 



„n— 3 



E 2 (x) 2(n-l) 



9 -x-\ 



2(n- 1) 



(n- l)(n + 2) 



n+2 
X 2 
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the equation for w has the following canonical form 

ww'(z) + w = <&(z) 

where 



— n — 6 n — 6 re — 2 re — 2 

i?((n - l)(n + 2))»+2 z«+2 - p ((n - l)(n + 2))™+2 z™+2 . 



The transformation formulas between (x, y) and (z, w) are given by 
(2.1) 



I n + 2 

X 2 



(ra-l)(ra + 2) 

re — 2 

id = yx 2 , 



and 

) x = ((n - l)(n + 2)z)^ 



(2.2) 

We summarize the above discussion as 



_ re— 2 

y = w ((n - l)(n + 2)z) n + 2 



Proposition 2.1. A gradient Yamabe soliton metric g on M is determined uniquely by one of the 
fallowings: 



(1) A positive solutions ip = <p(r) to equation (1.2). 

(2) A solution w = w(z) to the following equation defined on positive z-axis: 
(2.3) ww'(z) + w = ${z) 

with 

— re — 6 n—2 

$(z) =R((n- l)(n + 2)z)™+ 2 - p ((«- \){ n + 2)z)«+2 . 

(3) A trajectory a = (z(s),w(s)) on the half plane with positive z of the following dynamical 
system: 



(2.4) 



z{s) = w 
w(s) = <&(z) — w, 



where $(z) is given in (2) above. 

Moreover the metric g is complete if and only if either <p = <p{r) is defined for all r E K, or the 
following integral 

f 1 2 

(2.5) 1=1 -z'^dz 



w 

diverge in both directions along a trajectory a in (3). 

Proof. We already showed the equivalence of (1) and (2). The dynamical system formulation of 

(2) gives us (3). Next we consider the completeness of the metric g defined by ip or a trajectory a 

(3) . The metric g is complete if and only if its restriction on the real line R which is a geodesic in 
M is complete, i.e., ip{r) exists for all r G (—00,00) and stays positive. Suppose a = (z(s),w(s)) 
is defined for s £ (sq,si) where sq or s\ may be infinity and let (z^iUj) = ]hn s -> Si (z(s),w(s)) for 
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i = 0,1. Take another point (a,b) on a between (zo,wq) and (z±,Wi). Since dr = -dx which is 

a constant multiple of w z n+2 dz, the completeness of g is equivalent to the divergence of the 
integral I in both directions along a. Note that it diverges to infinity with different signs. □ 

Remark 2.2. In the rest of the paper, if two integrals I\ and 1% diverge or converge simultaneously 

along a trajectory, then we write I\ ~ I2. 

1 1 1 _ 1 2_ 

Remark 2.3. In the dynamical system (2.4) since ds = — dz and dr = w z ™+ 2 dz, we have 



(2.6) ds = z"+ 2 dr, 

i.e., s points in the positive direction of r and it is comparable with r when the trajectory approaches 
to a point [zq,wq) with zq > 0. 

Remark 2.4. If we use the variables x and y instead of z and u;, then the dynamical system is given 
by 

x(t) = 2(n — l)xy 

y(t) = R-(n-l)(n-2)y 2 -x\y + p), 
and the integral for the completeness of g is 

1 



(2.7) 



dx 

y 



for a trajectory a in xy-plane with positive x. 



Remark 2.5. The system (2.7) is more convenient in the case when cp has a zero point, i.e., the 
metric is rotationally symmetric. We give a brief argument showing that there is a unique complete 
rotationally symmetric metric g in this case based on the phase plane analysis. The initial condition 
that <p(Q) = and f'(0) = 1 indicates that we should look at the trajectories that approach to the 
point (0, 1) in the xy-plane. This point is a critical point as x = y = 0. The linearization of the 
system at this point is given by 

'x\ _ /2(n - 1) \ /x N 

jj \ _ 2 (n-l)(n-2)J [y y 

It follows that (0, 1) is a saddle point. One can see that two separatrices are the y-axis, and the 
separatrix 7 with horizontal tangent and positive x defines the Riemannian metric g. The trajectory 
7 approaches to another critical point as t — > 00 if p > and to infinity if p < 0. In either case, 
one can show that the integral J y~ l dx diverges as t — > 00 which gives the completeness of g. 



3. Gradient Steady Yamabe solitons 

In this section we prove Theorem i.e., any complete gradient steady Yamabe soliton on 
a warped product is necessarily a Riemannian product with constant <p. From the theorem by 
Cao-Sun-Zhang and Remark 1.5 in the introduction, we only have to show 

Theorem 3.1. Let n > 3 and (iV"^ 1 , y) be a Riemannian manifold with positive constant scalar 
curvature. Then there is no complete gradient steady Yamabe soliton metric on M n =lx v iV™ -1 . 



GRADIENT YAMABE SOLITONS 



7 



Our proof of the theorem above is separated into two cases, the case with dimM = 6, see 



Theorem 3.3, and the case with other dimensions, see Theorem 3.4 



First we fix the notions and state some general facts for all dimensions. Since M is steady soliton, 
we have p = and thus 

n — 6 n — 6 

3>(z) = where A = ((n - l)(n + 2))»+2 R > 0. 



By re-scaling the metric g if necessary we may assume that A = 1. So equation (2.3) and the 
dynamical system (2.4) have the following form: 

(3.1) w(z)w (z) + w(z) = z n + 2 , 

and 



(3.2) 
Let 




w 

n — 6 
Z n + 2 



w. 



V = {(z,w) G R 2 : z > 0} 



denote the half plane where we are looking for trajectories. We define the following curve S: 

S = {(z,w) G V : w = . 

It separates the half plane into two pieces: 

Vi = {{z, w) eV :w> , V 2 = {{z, w) G V : w < *(«)} . 

In any dimension we show that a trajectory does not define a complete Riemannian metric when 
it approaches the w-axis, the boundary of V. 



Lemma 3.2. Suppose a is a trajectory of the system (3.2) that approaches a point po = (0,wq) on 
the w-axis. Then the integral I along a does not diverge near p^. 

Proof. If po is not the origin, i.e., wq ^ 0, then in the integrand of /, \w~ 1 1 is bounded from infinity 

near the point po. Since ^2 < 1) tlie- integral I converges. If po is the origin, in the case when 

n = 6 we can use the explicit solution to the equation (3.1) as $(z) = 1, see for example |PZj . to 

see that I converges near the origin. In the follow ing we give a general argument for all dimensions. 
i i ^ 

Let u = z n + 2 as z > and then the equation (3.1 ) of z and w is equivalent to 

dw n + 2 
du w 

and the integral / is given by 



lu 



2n-5 



U 



n+l 



W 



I 



(n + 2) J 



ii 



n-l 



-du. 



The above equation of u and w defines the following dynamical system 



(3.3) 



u(t) 
w(t) 



(n + 2) (■ 



u 



2n-5 



U 



n+l 



W 



s 
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If n = 3, then the linearization of system (3.3) has the following form 



u(t)\ = (0 1\ (u 
w(t)J \5 OJ \w 

The above coefficient matrix has eigenvalues ±\/5 with eigenvectors (±1, v5). So the origin is the 
saddle point and there two different trajectories approach it, one with t — > — oo and the other with 
t — > oo. We use the power series to approximate these two trajectories. Let 

w = u k (ao + a\u + • • • ) 

with k > and ao 7^ 0. Then comparing the both sides of the equation on u and w shows that 
k = 1 and (Zq = 5. So for some a > small enough, we can estimate the integral I as 



(n + 2) f 
j 



— du ~ / — du = / udu < 00. 
w Jo u J 



If n > 4 then the linearization of the system (3.3) has the coefficient matrix 



1 


which shows that (0, 0) is a non-hyperbolic critical point of the nonlinear system. However since 
the coefficient matrix is not a zero matrix, by \Pe\ pp. 180, Theorem 2] we know that (0, 0) is a 
topological saddle, i.e., there are two different trajectories in V approaching this point. A power 
series approximation as in the case when n = 3 shows that w = u n ~ 2 (ao + a\u + • • • ) with nonzero 
clq around the origin and then / is finite when z — > 0. □ 

3.1. The case with dimM = 6. We have &(z) = 1 and S is horizontal line w = 1. The following 
Figure [T] shows the phase portrait of the dynamical system (3.2). Note that the horizontal axis is 



the z-axis and the vertical one is the u-axis. The green line is S and the half plane V is separated 
into two regions. 



Note that S is a trajectory of the system (3.2) and then any trajectory passing through a point 
in T^i (or V2) will stay in that region. 

Theorem 3.3. There is no complete steady Yamabe soliton on M e = M. iV 5 with N 5 being 
positive scalar curvature. 

Proof. To show the statement we only have to show that any trajectory cannot define a complete 
metric. Note that there is no vertical asymptote in the half plane V. Since along a trajectory we 
have w'(z) = ~ — 1, the u-axis is not a vertical asymptote of the trajectory. First the trajectory 
w = 1 does not define a complete metric. Next we consider the trajectory that lies in the region 
V\. Suppose 7 is a trajectory passing through a point 7(0) = (zo,uo) € V\. Then w(z) is a 
decreasing function in z and it meets u-axis with positive value of w. This shows that the metric 
is not complete as z — > 0. 

Now suppose (zo,u>o) 6 Vi and wq > 0. Since w is an increasing function in z, 7 meets either 



u-axis or z-axis as z decreases. If it meets w- axis (including the origin), then Lemma 3.2 shows that 
the metric is not complete when z — > 0. If 7 meets z-axis with positive value, then it intersects 
with z-axis vertically and thus enters the 4th quadrant. In this quadrant w is a decreasing function 
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FIGURE 1. Phase portrait with n — 6, A = 1 and p — 

in 2. Since u>'(z) = -Ep < —1 we have w(z) < —z for z large enough, i.e., — w > z. For the integral 
/ when z — > oo we have 

00 i _i r 00 i _i r 00 _s 

—z idz > / — z 4 dz = — z ^dz > — oo. 



So the trajectory does not give a complete metric. This finishes the proof. □ 

n — 6 

3.2. The other case with dimM 7^ 6. We have &(z) = z n + 2 and it is zero or approaches 00 as 
z — > when n > 6 or n < 6. The typical phase portraits in this case are shown in Figure [2] for 
n = 5 (the one on the left) and n = 8(the one on the right). The green curve is S for each case. 

Theorem 3.4. For any n / 6 there is no complete steady Yamabe soliton on M n =Ix„ N 71 ^ 1 
with iV" -1 being positive scalar curvature. 

Proof. The proof is similar to the one when n = 6. Since S is not a trajectory, any trajectory 
starts from a point on S will enter the region V\ or 7-V Note that there is no vertical asymptote 
of trajectories in the half plane V . When n > 6, the w-axis is not a vertical asymptote either. We 
claim that the w-axis is not a vertical asymptote when n < 6. Suppose not, there is a trajectory 
7 = (z,w(z)) in the region V\ such that along 7 we have 

lim w = 00. 

z^O 

Since in the region V2 with positive iu, w decreases as z tends to zero, 7 is bounded from below by 
the curve S. In particular it follows that there exists a sequence {zi > 0} that converges to zero, 
such that $'(zi) > w'(zi), i.e., 

re — 6 

— z< n+2 > - 1. 

n + 2 w[Zi) 
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FIGURE 2. Phase portraits with n = 5(the left one) and n = 8(the right one), A = 1 and p — 



Multiplying z n+2 on both sides yields 



n 



6 



> 



1 



, n+2 



n + 2 w(zi) 1 

Let i — y oo, then the right hand side converges to zero. However the left hand side is a fixed 
negative number which shows a contradiction. 

As in the proof of case n = 6, from the monotonicity properties of the function w = w(z) in 
different regions, one can show that a trajectory either meets the w-axis with finite w, or enters 
the region V2 with negative w and then decreases at least like the function w = —z for z large. 
In either case, the integral I is finite and so the metric defined by the trajectory is not complete 
which finishes the proof. □ 



4. Gradient Yamabe Shrinkers 



In this section we prove Theorem 1.3, the existence of complete gradient Yamabe shrinkers on 
warped products which are not isometric to the Riemannian product. We consider the case n = 6 



first and prove Theorem 4.8, and then the case with other dimensions and prove Theorem 4.13 



First we fix notions for this section and show some general facts for all dimensions. Let R > 

_, n — 6 — 

be the scalar curvature of N n ~ , A = ((n — l)(n + 2)) n + 2 and A = RA > 0. We may assume that 
p = ((n — l)(n + 2)) "+ 2 > by re-scaling if necessary and then we have 

n — 6 n—2 



So the differential equation (2.3) in z and w is given by 

. n— 6 n—2 

(4.1) w(z)w(z) + w(z) = \z™+2-z™+2. 
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and the dynamical system (2.4) has the following form 

J z(s) = w 



(4.2) 



n— 6 re — 2 

w{s) = Xz n + 2 — z n + 2 — w. 



First we observe that for all A and n, (£, 0) is a critical point of the dynamical system (4.2 ) where 

n + 2 

£ = A < . 

If n > 7, then (0, 0) is another critical point of the system. In the case of the steady solitons, we 
showed that the integral for completeness test is finite along a trajectory when it approaches the 
origin. Similarly we have 

Lemma 4.1. Suppose 7 is a trajectory in V that approaches the origin (0,0), then the integral I 
is finite. 



Proof. We use the variable u = z n + 2 around the origin (0,0). The differential equation (4.1) of z 
and w can be written as 

— = ( n + 2) (n + 2)u n+ . 

du w 



Following the argument of Lemma 3.2 in the steady case, we know that (0, 0) is a topological saddle 



of the dynamical system in u and w. A formal power series expansion shows that 

w = u n ~ 2 (ao + a\u + • • • ) 

where (n — 2)oq = (n + 2)A. So for some a > small we have 

f a 1 f a 1 

I = / — dz ~ / „ dz < 00 

Jo wz^+z Jo Z"+ 2 

which finishes the proof. □ 

The next two propositions characterize the local behavior of the trajectories around the critical 
point (£, 0). 



Proposition 4.2. The nonlinear system (4-2) has (£, 0) as either a stable node when (n + 2)A > 16, 
or a sable focus when (n+2)A < 16. In both cases the integral I diverges to infinite along a trajectory 
that approaches this point. 



Proof. At the point (£,0), the nonlinear system (4.2) has the following linearization 

(i) - {-X -J (; 

Let 7*1 and r2 be eigenvalues of the coefficient matrix of the above linear system. We have 

4 

r\ + r 2 = -1 and r x r 2 = -, — -^tv- 

(n + 2)A 

If (n + 2) A — 16 > then the two eigenvalues are negative real numbers and so (£, 0) is a stable node 
of the linear system. If (n + 2)A — 16 < 0, then the two eigenvalues are complex with negative real 



part and so (£,0) is a stable focus. On the right hand side of the nonlinear system (4.2), we only 
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have power functions in z and w with positive exponents, from |Pe[ pp. 142, Theorem 4] the critical 
point (£, 0) is a stable node or a stable focus of the nonlinear system respectively. In both cases if 
7 is a trajectory around this point, then 7 approach (£, 0) as s — > 00. Since the distance function r 
is comparable with the parameter s near the critical point, r tends to infinite as 7 approaches this 
point. □ 

In general we only know that a trajectory near the point (£, 0) will approach it as s — > 00. In 
this case we can show the convergence for a large area. 

Proposition 4.3. For any < zq < £, the trajectory with initial point (zo,0) converges to the 
critical point (£,0) as s —> 00. 

Proof. Suppose 7 is a such trajectory with 7(0) = (zq,0). As s becomes positive, 7 enters the 
first quadrant and then w = w(z) is an increasing function of z. As z(s) = uu > 0, z and uu are 
increasing as s increases. Then 7 will meet Si and after that z increases while w decreases. The 
trajectory 7 will meet the z-axis and enter the 4th quadrant. Then it will meet Si and the z-axis 
again. Suppose 7 meets the z-axis from the 4th quadrant at s = si, i.e., 7(si) = (zb,0). 

We claim that Zb > zq. Suppose not, since there is no periodic orbit by Bendixson's criteria, see 
for example [Pet pp. 245, Theorem 1], we have Zb < zq. Then the curve 7(5) with s G [0, si] and 
the segment {(z, 0) : Zb < z < zq\ form a bounded region and 7(5) stays in this region when s < 0. 
So 7(5) will approach a critical point as s — > —00. This contradicts the fact that in this bounded 
region there is only one critical point (£, 0) which is stable. 

Now we have Zb > zq and we have a bounded region E whose boundary consists of the curve 
7(5) with s G [0, si] and the segment {{z, 0) : zq < z < Zb}- The trajectory 7(5) for s > si stays in 
E and so 7(5) approaches the unique critical point (£,0) as s — > 00. □ 

In the following we show that there is a unique trajectory 7(5) = (z(s),w(s)) that approaches 
(£, 0) as s — > 00 and the function w = w(z) has a distinguished asymptotic expansion at infinity 
when s — > —00. Such asymptotic expansion of w implies that the integral I diverges when s — > —00 
and so the metric defined by 7 is complete. 

4.1. The case with dimM = 6. We have 

= A - y/z, 

and 

S 1 = {(z,w) eV:w = X — \[z\ . 

The curve S\ determines the monotonicity of trajectories. We introduce another curve in the half 
plane V which determines the convexity of trajectories. From the formula w"(z) of a trajectory 
a = (z,w(z)) we define 

S 2 = l(z, w) € V : ^ - 2(A - ^fz)w + 2(A - ^fzf = J . 

The curve Si separates the phase plane V into two pieces: 

y>\ = {(z,w) € V : w > A - y/z] , V 2 = {{z,w) G V : w < A - xfz] . 



GRADIENT YAMABE SOLITONS 



13 



Note that Si is not a trajectory of the system in z, w. Any trajectory starts from Si will enter one 
of these two regions. 
Let 

fx = (A - \fz)\[z + (y/z - \)\J z - 2yfz 

h = (A - yfz^yfz - {yfz — \)\Jz~- 2yfz 

for z > max {4, A 2 }. On the 4th quadrant, S2 has two components given by different formulas: 

S2a = {(z,w) £ V : w = fi, z > max(4, A 2 )} 

S 26 = {(z,w) e V : w = f 2 , z > max(4,A 2 )} . 

These two curves have a common point T when z = max {4, A 2 }. If A > 0, T is the same as the 
critical point (£, 0). If A < 2, then T has the coordinate (4, 2A — 4). It is easy to see that when z 
is large, the defining function of S2 a is asymptotic to w = —yfz and the one for S26 is asymptotic 
to w = —2z. We define the following trapping region T in Vi'- 

T={{z,w)eP 2 :f2<w<f 1 }. 

The following Figure [3] shows the region T for some typical values of A. The green curve is Si 
and the dotted blue curve is S2. 




FIGURE 3. The region T for A = 1,2 and 3 when n = 6. 

Proposition 4.4. For these two curves Si and S 2 we have 

(1) If A < 2, then on S2 we have z > 4 and the whole curve S2 is in the J^th quadrant as 
S2 = S2a U S26 and it is below Si . 

(2) If X > 2 then S 2 has non-empty components in the 1st quadrant for 4 < z < A 2 , and in the 
4th quadrant for z > A 2 which is S 2a U S2b- The curve S2 is above Si in the 1st quadrant, 
and is below Si in the 4th quadrant. 

Furthermore if a trajectory meets the boundary ofT with intersection point (zq, wq), then it enters 
and stays inside the region T for z > zq. 
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Proof. The defining equation of S2 is given by w"(z) = and we can rewrite it as 
w"(z) = — ((«, - (A - v^)v^) 2 - v^(v^ - 2) (A - V~zf) ■ 



It is easy to see that w = implies that z = A 2 for z > 0. So the set of = with positive z 

consists of the single point (A 2 ,0) and the curve S2. At any point on S2 we have 

J~z(J~z - 2)(A - ^/z") 2 = (iu — (A — \fz)\fz) 2 > 0. 

It follows that -y/z — 2 > 0, i.e., z > 4. Moreover the solution to the equation = is given by 

wi,2 = (A — \fz)\fz ± |a/z — A| \j~z~^-2\fz. 
If A < 2, then ^/i" - A > and 

wi,2 = - A) {-^fz ± ^z-2^) < 0, 
i.e., S2 lies in the 4th quadrant. If A > 2, then we have 

2 = (A - Vi) (^/z ± yj z - 2y/z^ > if4<z<A 2 , 



and 



u>i,2 = (v^ - A) ^— \/i ±\Jz- 2y/zj < if^>A 2 . 

So the curve S2 has nonempty components in the 1st and 4th quadrants. Using the formulas of 
wi t 2 the relative position between Si and S2 can be checked by the sign of A — \fz — w\p when 
z > A 2 (and 4 < z < A 2 if A > 2). 

Next suppose a is a trajectory that touch S^a & t the point p = (z, w). Since = w < in the 
4th quadrant, z increases as s decreases. Let mi = w'(z) be the slope of the tangent of a at p and 
^2 = f[( z ) be the slope for the curve S^a- We have 

A — Jz A - Jz 1 1 

mi = — 1 = 7 r- - 1 = - 1 

h (A - Vi) [Vz - \A ~ 2 Vz) ^ ~ V z ~ 2 v^ 

2^2^ 

m 2 = -1 + — = + V - v + ^ l _ 

2^i 2^ 2y/z-2yfi V 



and then 



A 1 y/z - A / 1 A ^ - A / - 1 
m 2 -mi = — -7= — - H . = 1 



2^ 2 2y/z-2y/z V v^/ 2^/i \^jz-2y[z~ 
> 0. 
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So a will enter the region T as z increases and it cannot escape from this region through the curve 
as z increases. Suppose a touches the curve at p = (z, w). Let m,i and m-2 be the slopes of 
the tangents of the curves a and S2& at p. Then similarly we have 

1 ^z - 2jz 

2 2^/z 



and 



A y/z -2yfi v^-A / 1 

mo = — 1 H 7= — 7= , 1 — 

2^z 2jz 2 ^z -2^\ Vz 



mi - m 2 = \ r- / +1 > 0. 



\y/ Z -2yfi 

So a will enter the region T as z increases and it cannot escape from it through the curve 52b- 
It follows that if a meets the boundary of T, then it will enter and stay inside this region as 

z — > CO. □ 

We would like to study how the trajectories behave at the infinity. Let u = y/z, then the equation 
of z and w is equivalent to the following one in u and w: 

(4.4) — = — ^ 

au w 

The above equation defines the following dynamical system 

u(t) = w 

w(t) = 2u(\ — u — w). 

It has (A, 0) as a critical point which is a stable focus(if A < 2) or a stable node(if A > 2). The 
system of u, w has one more critical point at (0,0). The linearization shows that it is a saddle 
point. 

The rough picture of the behavior of the system at infinity can be seen from the phase portrait 
on the Poincare sphere by the polar blow-up technique, see for example |ALGMj . Namely, we 
introduce the new coordinates X, Y and Z such that 

X Y 

u= z> w= z- 

Then the phase portrait of the new system in X, Y, Z near the equator X 2 + Y 2 = 1 of the Poincare 
sphere § 2 = {(X, Y, Z) £ M 3 : X 2 + Y 2 + Z 2 = 1} characterizes the behavior of the origin system 
in u, w at infinity. One can see that (1,-1,0) is a critical point and it is a saddle-node. The 
separatrix between the two hyperbolic sectors approaches the focus/node point (A, 0) (if A > ^) or 
the saddle point (0, 0) (if A = ^) in ii-u;-plane. If A < i, it approaches other critical points (0, ±1, 0) 
on the equator of the Poincare sphere. When A > \ this separatrix corresponds to the unique 
trajectory 7 in uw-plane such that along 7 we have 

lim (u(t),w(t)) = (A,0) or (0,0) lim u(t) > and lim ^| = -1. 

t— >oo t— >— 00 t— >— 00 wit) 

Note that when A = \, 7 is the graph of the function w = —u or w = —\fz. 
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In the case when A > ^, if the trajectory 7 has positive u for every t, then there is a corresponding 
trajectory 7 in the zw-plane with positive z approaching (£, 0) as s — > 00 and hin^^oo -7= = — 1. 
One can see that the integral / is unbounded as z — > 00, i.e., the metric given by the trajectory 7 
is complete. 

The above rough picture indicates that we can show the existence of a complete Yamabe shrinker 
for A > \ in two steps: 

(1) There exists a trajectory 7 along which, the function w = w(z) is asymptotic to w = —\fz 
for z large(or s — > —00). 

(2) When s — > 00, 7 meets the z-axis with < z < £, and then approaches (£,0). 

We verify the statements above in the following two lemmas. 



Lemma 4.5. There exists a positive number Zq and a trajectory 7 of the system (4-2) in z,w such 
that 

(1) 7 is in the 4th quadrant when z > Zq, and 

(2) along 7 we have 

lim —= = — 1. 



Proof. Let x = ~ and y = ~, then the equation (4.4) becomes 

(4.5) x 2 ^- = 2Xxy 3 - 2y 3 - 2y 2 - xy 

ax 

We consider a formal power series solution to the above equation 

y = x k (a + a\x + a 2 x 2 H ) 

such that ao < 0. Comparing the lowest term on both sides yields two cases. In the first case k = 1 
and the recurrence relation of ai is given by 

ao = —1, a\ = —2, A = 2 

and 

(i + 2)oj = 2A 

a p a q a r — 2 ^ a p a q — 2 ^ a p a q a r , for % > 3. 

p+q+r=i—2 p+q=i p+q+r=i—l 

So this case can only exist when A = 2 and then there is a family of formal power series solutions 
parameterized by 02 £ M. 

In the second case we have k = and the recurrence relation of aj is given by 

a = -1 

and 

(z + l)dj = 2A a p a q a r — 2 a p a q — 2 a p a q a r , for i > 0. 

p+5+r=i p+q=i+l p+q+r=i+l 

Solving aj+i from the above relation yields 

2aj+i = — (i + l)oj + 2A a p a q a r — 2 — 2 a p a q a r , 

p+q+r=i p+q=i+l p+q+r=i+l 



GRADIENT YAMABE SOLITONS 



17 



where the indices p, q and r in the last two sums are positive integers. It follows that all the 
coefficients 04 are determined, i.e., the formal power series in this case is unique. We apply [Wa, 
Theorem 33.1] to the equation (4.5) and we conclude that there exists an e > such that for all x 
with \x\ < e there is at least one solution y = 4>(x) admits the asymptotic expansion 



4>{x) ~ — 1 + o,{X l , x — > 0. 

i=l 

If we restrict x to be positive and translate back to the variables z and w, then it says that for 
large z > 0, there is a solution w = w(z) to the equation (4.1) that is asymptotic to the function 
w = —y/z. So we can choose Zq > such that w = w(z) < for any z > Zq and the trajectory 7 
is the one that passes a point in the graph of the solution. □ 

Remark 4.6. The technique of asymptotic expansion for a differential equation is also used by R. 
Bryant to show the existence of complete rotationally symmetric steady Ricci soliton metric on M. n , 



sec Bry|. 

Lemma 4.7. Suppose 7 is the trajectory in the previous lemma. Then 

(1) it stays inside the half plane V , 

(2) it stays between Si and S2a in the J^th quadrant, and 

(3) it approaches the critical point (£, 0) as s — )• 00. 
Moreover such trajectory with the above properties is unique. 

Proof. We already knew that if a trajectory passes a point below the z-axis and above S±, then 
it will meet the z-axis and then the value of w stays positive. Similarly since w"(z) < in the 
region T C V2, any trajectory passes a point in T will stay in this area as z increases and then w 
decays at least as fast as the linear function w = — z as z — > 00. Since along 7 the function w is 
asymptotically like the function — y/z, the trajectory 7 stays between Si and 62a- 

If A > 2, then Si and meet at (£ = A 2 , 0). Since z = w < 0, z decreases as s increases. There 
is no critical point other than (£,0) in the region bounded by Si and S2 a , so 7 will approach this 
point as s — > 00. 

Now we assume that | < A < 2. We consider another curve S3 in the 4th quadrant: 

S 3 = {( z ,h(z)) : h(z) = -y/z with z > 0} . 



Suppose a is a trajectory of the system (4.2) and it touches S3 at (z,w). Then the slope of a at 
the intersection is given by 

A-v^ , A 
m = -=- - 1 = — T 

-y/Z yJZ 

and thus we have 

It follows that if a lies below the curve S3 for some zq > 0, then it stays below S3 for all z > zq. 
Next we consider the relative position between S3 and S2a- Since 



fl( Z ) - h( Z ) = (A - y/z) fy/z - \Jz-2y/z^j + ^fz ^ X - ^ 



> 0, as z — > 00, 
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it follows that lies above S3 for z > large. In particular if a trajectory lies below S3 for some 
zq > 0, then it will stay below the curve S2 a for ah large z. So the trajectory 7 stays in the region 
E bounded by Si, S2 Q , S3 and the 2-axis. We consider the convergence of 7 when s tends to 00. 
It will meet the z-axis at (zq, 0) with zq > 0. We claim that zq > 0. If not, then 7 is a trajectory 
emanating from the origin. Suppose on 7 we have w = w(z) for small positive values of z. From 



the differential equation (4.1) we have the following series expansion 

w(z) = — v2X^/z + higher order terms. 
Since A > i, for small values of z the trajectory 7 stays outside the region E which contradicts 



the fact that 7 lies inside E. From Proposition 4.3 we conclude that 7 approaches the stable focus 
(£, 0) within the half plane V as s — > 00. 

To finish the proof, we show that such trajectory is unique. Suppose not, then there are two 
solutions 101(2) and 102(2) which define two trajectories 71 and 72 with the stated properties in this 
lemma. Since they cannot cross each other in the 4th quadrant, we may assume that < 102(2) < 



101(2) for all z > zq with some positive Zq. From the differential equation (4.1) we have 



, , . A — yfz A — sfz A — yfz , . 

W 1 (2) - W 2 (Z) = = (Wl - 10 2 ) • 

Wl W2 W1W2 

Since both 71 and 72 stay between Si and S2 a , the functions 101(2) and 102(2) are asymptotic to 
the function —\fz for large z. However from the above differential equation and the comparison of 
first order differential equations, i.e., the Gronwall's lemma, we see that log(ioi — W2) is asymptotic 
to the function 2y/z which shows a contradiction. □ 

From the previous results we obtain 

Theorem 4.8. Let M 6 = R x„ iV 5 and N has constant scalar curvature R > py/W. Then there 
exists a unique complete Yamabe shrinking soliton metric on M with constant p that is not isometric 
to the Riemannian product. 



Proof. For any A > \ with p = 2 ^^ , Lemma 4.5 and Lemma 4.7 imply that there exists a trajectory 



7 of the dynamical system (4.2) such that 

(1) z is positive along 7, 

(2) 7 approaches (£, 0) with £ = A 2 as s — > 00, and 

(3) it satisfies the following asymptotic properties: 

lim z = 00, lim ^= = —1. 

Property (1) ensures that <p is positive everywhere, i.e., if defines a Riemannian metric g = dr 2 + 
<p 2 (r)gj\r on M 6 = lx ip iV 5 . Since the Riemannian distance r is comparable with s near the point 
(£, 0), the metric g is complete when s — > 00. Using the asymptotic properties of 7 we have 

roo poo g 

1= r dz ~ — / z~±dz = — 00, 

Ja WZ4, Ja 

for any a > 0, i.e., the metric g is complete as s — > —00. So 7 defines a complete Yamabe shrinker 
on M. Under the re-scaling of the metric g, R/p = X/p is unchanged. So we have the scalar 
curvature R of iV 5 bigger than p\/W. □ 
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Remark 4.9. From the proof in Lemma 4.7 one sees that the trajectory 7 is bounded by the curve 



w = —\fz from below. It implies that the scalar curvature R of the metric on M defined by 7 is 
positive. To see this, let p = ^= and then we have 



„ / w 1 / w v 

R = p + v / = p + y = p + —— = l + > 0. 

y/AQyfz V40 V yjz; 

4.2. The other case with dimM 7^ 6. In this case we have 

n — 6 n — 2 

$(z) = Xz n + 2 - z«+2 . 



The two curves that will be useful to study the trajectories of the dynamical system (4.2) in z and 
w are 

5i = {(z,w) G V : w = $0)} 

and 

$ 2 (z) 



The half plane 7-" is separated by S\ into two parts: 

Vi = {(z, w) eV :w> *(z)} , = {(«, tu) G V ■ w < *(«)} . 

When n < 6 let z a be the unique positive solution to the equation 4<5'(z) + 1 = 0, or equivalently 

n-2 ^ n + 2 4, 
n — 6 4(n — o) 

When n > 6 and (n + 2)(n — 6) A < (n — 2) 2 , the above equation has two distinct positive solutions 
and let z a be the larger one. If n > 6 and (n + 2)(n — 6)A > (n — 2) 2 , then the above equation has 
no real solution and we let z a = 0. It is easy to see that in this case (£, 0) is a stable node. In these 
three cases, we define two functions for z > max{z a ,^}: 

We will see that both functions are non-positive and fi(z) > / 2 (z) for any z > max{z a ,^}. We 
call the graphs of fi^ the curves S^a and 526, i«e., 

5*2a = { (z, w) G P : w = /1 (z) , z > max(z a , ()} 

= {(z,w)eP:w = /2(4 z > max ( z «.0}' 

Proposition 4.10. Lei n < 6. For i/iese too curves Si and S2 we have 

(1) If (n + 2)A < 16, i.e., (£, 0) is a stable focus, then on S2 we have z > z a and the whole 
curve S2 is in the J^th quadrant with S2 = S^a U ^26 and it is below Si. 

(2) // (n + 2)A > 16, i.e.,(£,0) is a stable node, then S2 has non-empty components in the 1st 
quadrant with z a < z < £ and in the J^th quadrant with z > £. In the 1st quadrant S2 is 
above Si, and in the J^th quadrant 5 2 = S^aUS^ft an d ^ s below Si. Moreover the three curves 
Si, S2a and 5*26 have the unique common point (£,0) 
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Proof. The defining equation of S2 is given by w"(z) = and we have 



w 2 V J w 2 \\ 2<S>'(z)J A($>(z)) 2 ' 

So on the curve S2 we have 1 + A&{z) > 0. We consider the monotonicity of &'(z). Since 

1 8 / 4 \ 

&(z)= +2 -z "+2 (^(n - 6)A - (n - 2)z«+2 J 

and 

= ^ n + 2 y z 11+2 ( 4 ( n - 2)^"+ 2 - 8(n - 6)AJ > for z > 0, 

the function <5'(z) is an increasing negative function from — \ and bounded above by zero as it 
approaches zero when z — > 00. It follows that yjl + 4$'(z) is an increasing non-negative function 
which is bounded above by 1. Let z a > be the unique value of z such that <£>'(z a ) = — \ and then 
we have 

n - 2 _l_ n + 2 _^^. r 

a Zn+2 ~ 77 a\ Zn+2 - A for Z - Z <*> 

n — b 4(n — b) 

with equality if z = z a . So the defining equations of S2 can be written as 

and the domain of u+2 is [z a , 00). If u>i or u>2 has positive value for some z > z a , then &(z) > 0. 
It follows that z a < £ which is equivalent to the following inequality 

n - 2 n + 2 x2 

« A " 77 S7 A > A ' 

n — b 4(n — b) 

i.e., (n + 2)A > lb. 

Next we consider the relative position between S\ and S2. Note that when z > max{z a ,^} we 
have wi ; 2 = fi,2- If (w + 2)A < lb, then S2 is in the 4th quadrant and fi > f2- Since z a > £, 
$(z) < for z > z a and then we have 

h 1 /^77^\ 2 



*(*)" 2$'(z)( 1 V / T+4$ 7 M) - 1 + v/l + 4$/(z) 

So the curve 5i is above 52 in the 4th quadrant. If (n + 2)A = lb, then a similar argument shows 
that the two curves S\ and S2 have a unique common point (£, 0). If (n + 2) A > lb then we already 
know that Si is above the curve 52 in the 4th quadrant. We consider the curves in the 1st quadrant 
where z a < z < £. Since <&(z) > 0, we have 

|i = (1 ± yiTi^) = T - 7 J_ > 1, 

i.e., 52 is above the curve Si. Note that in this case Si and 52 also have a unique common point 

(£,0). □ 



> 1. 
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In the case when n > 6, if (n+2)(n — 6)A < (n — 2) 2 then S2 has two connected components in the 
half plane V . The one given by ^U^j lies entirely in the 4th quadrant. If (ra+2)(n— 6)A > (n— 2) 2 , 
then 5*2 has only one component and intersects with Si at (£, 0). Furthermore we have 

Proposition 4.11. Let n > 6. For these two curves S\ and S2 we have 

(1) If (n + 2)A < 16, i.e., (£, 0) is a stable focus, then the restriction of S2 to {z > max(z a ,( t )} 
is S2a U S26 which lies entirely in the J^th quadrant and is below S\. 

(2) 7/(n + 2)A > 16, i.e.,(£, 0) is a stable node, then the restriction of S2 to {z > £} is S2 a US26 
which lies entirely in the J^th quadrant and is below Si. Moreover the three curves Si, S2 a 
and S2b have the unique common point (£, 0). 

Proof. The argument is similar to the one in the previous proposition. Note that when z > 
max {zq,, £}, we also have $"(z) > 0. □ 

From the formulas of /i,2 we see that fi = fi when z = max{z a ,£} and /i(or f2) is asymptotic 
to the function $(z)(or — ^zf^)- We define the following trapping area in V2' 

T = {(z, w) G V : f2 < w < fi and z > max(z a , £)} . 



Proposition 4.12. Suppose 7 is a trajectory of the system (4-2) with n ^ 6. // 7 meets the 
boundary of T, then it enters and stays in this area as z increases. 

Proof. Suppose 7 meets the boundary of T at (z, w). Let mi and m,2 be the slopes of the tangents 
of 7 and the boundary curve at this point. If 7 meets S2 a , then we have 

Hz) , Hz) 

w fi{z) 
and m,2 = f[(z). A straightforward computation shows that 

mi ~ m2 = J" { r^^ h ■ iM s 1 - (V 1 + - 1 + + m*) - 3) 

{$'(z)) 2 y/l + 4$'(z) a/1 + 4$'(z) -IV 
< 0. 

So the trajectory 7 cannot escape the area T through the curve S2 a - Similarly if 7 meets S2& then 
we have 



H Z ) , v 

mi-m 2 = y-^r- 1-/2(2) 

- f'jfi^ = (Vl + 4^(z) + l + ^)(Vl + 4^(z) + 3)) 

(^'(z)) 2 ^! + 4$ / (z) a/I + 4$'(z) + 1 v / 

> 

i.e., 7 cannot leave this area through the curve S2b either. So 7 will stay inside T as z increases. □ 

We have the following existence result of Yamabe shrinkers. 

Theorem 4.13. Let n ^ Q and p > 0. Let R > p^ ((n - l)(n + 2))^ an d (jV" -1 ,<7) be Rie- 
mannian manifold with constant scalar curvature R. Then there exists a unique complete Yamabe 
shrinking soliton metric on M n = K x^N 71 ^ 1 with constant p that is not isometric to the Riemann- 
ian product. 
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Proof. Recall that A = R((n — l)(n + 2))^+ 2 and p = ({n - l)(n + 2))~^+ 2 . The inequality of R 
and p in the theorem is equivalent to A > ^| by re-scaling. As in the case when n = 6 we show 
the existence and uniqueness of the metric in two steps. 

Step 1. We show that there is a solution w = w(z) for large z > with the asymptotic 
expansion: 

n— 2 

w{z) ~ — z n + 2 + lower order terms. 
Then we let 7 be the trajectory that passes a point in the grap h of this solution. 



_ 1 -1 

Using the variables x = z n + 2 and y = — , the equation (4.1) has the following form 



x 2n+i^ = ( n + 2 ) ((Ax 4 - l)y 3 - x n - 2 y 2 ) . 

For some e > small it has a solution y = 4>{x) with < x < £ with the following asymptotic 
expansion 



(x) ~ -x n ~ 2 + x n - 2 Y,a iX \ x^O. 



i=l 



It follows that the equation (4.1) admits a solution w = w(z) such that 

w(z) 



lim 



00 z^ 



Step 2. We show that the trajectory 7 obtained in the previous step, i.e., it is defined by the 
solution w = w(z), is unique, stays inside the half plane V and approaches to the critical point 
(&0). 



From Proposition 4.12 7 stays between the curve S\ and S^a- in the case when (£,0) is a 
node, i.e., (n + 2)A > 16, since S\ and S2 a meet at (£,0), 7 approaches this point as z decreases. 
The metric defined by 7 is complete at z = £ and the completeness as z — > 00 follows from the 
asymptotic property of the solution w = w(z). 

In the case when (£, 0) is a stable focus, we consider the auxiliary curve S3 in the 4th quadrant 

n— 2 

which is the graph of the function h(z) = —z n + 2 . One can show the following properties and the 
proof is similar to the case when n = 6: 

(1) If a trajectory a lies below £3 for some zq > 0, then it will stay below S3 for all z > zq. 

(2) The curve S*2 a stays below the curve S3 for all large z > 0. 

(3) If a trajectory emanates from the origin, then it will stay below S3. 

It follows that the trajectory 7 stays inside the region E bounded by the curves Si, S2 a , S3 and 



the z-axis, and it meets the z-axis with positive value of z. From Proposition 4.3 the trajectory 
7 stays inside the half plane V and approaches the stable focus (£,0). So 7 defines a complete 
Yamabe shrinking metric on M which is not isometric to the Riemannian product. The uniqueness 
of such metric, or equivalently of the trajectory 7, follows from a similar argument as in the n = 6 
case. □ 

Remark 4.14. The fact that the scalar curvature R of t he metric defined on M by the trajectory 7 



is positive can be seen from the proof in Theorem 4.13 Let p = ((n — l)(n + 2)) ™+ 2 and then we 
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)>„ 



Z n + 2 . 



n-2 
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